Critical Scattering and Dynamical Scaling in an Heisenberg Ferromagnet 
Neutron Spin Echo versus Renormalization Group Theory 
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High resolution Neutron Spin Echo (NSE) spectroscopy was used to investigate the dynamics of 
an 3D Heisenberg ferromagnet in the exchange-controlled regime over a broad range of temperatures 
and momentum transfer. These results allow for the first time an extensive comparison between the 
experimental dynamical critical behavior and the predictions of the Renormalization Group (RG) 
theory. The agreement is exhaustive and surprising as the RG theory accounts not only for the 
critical relaxation but also for the shape crossover towards an exponential diffusive relaxation when 
moving from the critical to the hydrodynamic regime above Tc- 
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Critical phenomena are ubiquitous in condensed mat- 
ter physics pQ and their theories are now used to under- 
stand the behavior of many complex systems, like sand- 
piles [2 or tropical rains [3]. Phase transitions are dy- 
namical processes in essence, specially the continuous 
phase transitions. These are characterized by diverging 
fluctuations at the critical temperature Tc, with a crit- 
ical slowing-down and a scale- free behavior: power law 
spatial and temporal evolutions of the order-parameter 
fluctuations. In the critical region, the dynamics is gov- 
erned by a single dimensionless variable, according to the 
dynamical scaling hypothesis [4] . Ferromagnetism serves 
as a paradigm for continuous phase transitions and as 
such a lot of theoretical and experimental effort [5j |6] 
has been devoted to understand and predict in details its 
static and dynamic critical behavior. 

Critical dynamics is a much more complex problem 
than the static critical phenomena as it depends not only 
on the static universality class but also on conservation 
laws and on the equations of motion of the system [7]. 
For example, the 3D Heisenberg ferromagnet displays 
two different behaviors, depending on the relevant in- 
teractions, exchange or dipolar. The balance between ex- 
change or dipolar interactions is quantified by the dipolar 

wave- vector Qjj (Q 2 D — Mo^b^jI - ' ^ 2 being the second 
moment of the exchange interactions). Exchange con- 
trols the dynamical behavior for Q > Qjj and the order- 
parameter (i.e. the magnetization) is conserved (model 
J of Reference ([7])). The critical dynamics at Tc is 
governed by a dynamical exponent z= D +^~ ri ; where D 
is the space dimension and 77 the Fisher exponent. For 



3-D Heisenberg systems, z « 5/2 (figure [lj. Dipolar in- 
teractions, on the other hand, control the dynamics at 
Q < Qd- These interactions do not conserve the order 
parameter and at Tc the dynamical exponent crosses over 
from z=5/2 (exchange) to z=2 (dipolar) at Q « Qjj/10. 
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FIG. 1: Critical slowing down of the characteristic relaxation 
time t c of the magnetic fluctuations in CdCriS^. The solid 
line is a power law fit, r c (Q, Tc) = ToQ~ z with z = 2.47±0.02 
and to = 0.02813 ns/nm~ z . 

In the critical region above Tc and in the sole presence 
of exchange interactions, the relevant scaling variable is 
x = k/Q, with n the inverse magnetic correlation length. 
In this region the magnetic fluctuations decay faster than 
an exponential with a shape cross-over to an exponential 
relaxation when going from the critical to the hydrody- 
namic regime for x ^> 1. This simple scaling hypothesis 



2 



breaks down in the presence of dipolar interactions. At 
Q < Qd, the relaxation time depends on two reduced 
variables: x = n/Q and y = Q/Qd [H] whereas the relax- 
ation recovers a simple exponential decay. The boundary 
conditions for the exchange regime are therefore fixed by 
the theory: Qjj < Q < ft, which implies that the criti- 
cal non-exponential relaxation can only be observed in a 
precise temperature and Q range and for ferromagnetic 
systems with low enough values of Qd- 

Independent renormalization Group (RG) [9] and 
Mode Coupling (MC)[5 calculations predict similar 
shapes for the spectral function at the exchange con- 
trolled regime at Q > Qd- These predictions, however, 
were never tested thoroughly. One reason ist that in neu- 
tron scattering experiments the finite resolution blurs out 
the magnetic relaxation spectra. Triple-axis neutron ex- 
periments (see [10] and references therein) focussed on 
constant-energy Q-scans which displayed a broad maxi- 
mum the position of which was compared with the the- 
oretical calculations. In this letter, we present the first 
detailed analysis of the shape of the relaxation function 
in a model 3D Heisenberg ferromagnet by taking advan- 
tage of the recent progress in the field of high resolution 
Neutron Spin-Echo (NSE) spectroscopy. 

The RG theory expresses the spectral shape function 
F{Q,uj) as the real part of a complex function of the scal- 
ing variable n/Q and the scaled frequency co/co c , with 
the help of an e = 6 — D expansion [9 and three pa- 
rameters /c, a (for the critical behavior) and b (for the 
cross-over to the hydrodynamic regime) : 



F{Q,u) = —9? 



(Z(x)Hi (x, iu)) 1 — i 



(i) 



The function TLi(x,iu) is the self-energy of the dynamic 
susceptibility , approximated by : 



III (x,iu) = ((l + bx 2 



(2) 



(where u = uo/uo c x Z(x)/X(x) for clarity). Z(x) is es- 
sentially a scaling function, approximated by: 



Z(x) 



(l + 6a; 2 )T 



1 — k x Arctan 



( a (l+6x 2 ) 2 ) 



(3) 



and X(x) is the x(Q) susceptibility modelled by a 
Lorentzian 1/(1 + x 2 ) in the Ornstein-Zernicke approx- 
imation. The scaling function describing the crossover 
with x from the critical behavior to the diffusive one in 
the hydrodynamic regime is given by : 



3 (Q,T) _ r c (Q,T c ) _ 1 Z(x) 



u c (Q,T c ) r c (Q,T) X(x)Z(0) 



(4) 



For a better comparison with the NSE experimental re- 
sults, a analytical form of the intermediate Fourier trans- 
form of the spectral function [IT] has been proposed: 



I RG (Q,t) =e- a(x) ^(Cos(b(x)-) +Sin{b(x)-)) (5) 
with 

2 93 -4- 3 9fi^» -4- 4 1 QF^ 2 9 Q3 

(6) 



, x 3.23 + 3.26x + 4.195x 2 x 



1 + 0.2SX 1 / 2 + 3.84x 3 / 2 



l + 2.1x 



and r c (Q,T) oc Q(x)tqQ~ z . We will use this conve- 
nient parametrization of the RG theory to reduce our 
NSE experimental data, at the critical temperature Tc 
(Sl(x) = 1) and above (Q(x) ^ 1). 

CdCr<iS± is a normal spinel, with all Cr 3+ ions hav- 
ing a magnetic moment of spin only type (S=§) with a 
very low magnetic anisotropy. CdCr^S^ is a semicon- 
ductor and undergoes a ferromagnetic transition below 
Tc = 84.8 K. The magnetic interactions are positive 
between the six nearest neighbors (Ji/fc^ = 13.25 K), 
negligible between the six next nearest neighbors and 
negative (J^/ks = —0.915 K) between the twelve third 
nearest neighbors. These interactions lead to a very 
small dipolar wave- vector of Qd = 0.51 nm" 1 , making 
this coumpond an ideal candidate for the study of the 
exchange governed critical dynamics. Our comprehen- 
sive investigation of the static critical properties ([12] 
and references therein), determined the critical expo- 
nents (3 = 0.33±0.03, v = 0.70±0.03 , rj = 0.065±0.030, 
and 7 = 1.40 ± 0.04, which are very close to the theoret- 
ical values for the 3D Heisenberg model [13]. Therefore, 
CdCr2S± is a textbook example of the 3D Heisenberg 
ferromagnetic class. 
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FIG. 2: Magnetic critical scattering at Tc in CV/CV2S4 , 
measured on the NSE spectrometer IN11 (ILL). Lines are 
adjustments of the RG theory of I(Q,t) (see equation [5| 

The measurements were performed on the neutron 
spin-echo (NSE) spectrometer IN11 at the Institut Laue 
Langevin( Grenoble, France). The polycrystalline sam- 
ples were enclosed in a sealed aluminium can and placed 
in a standard Orange cryostat, where the temperature 
was controlled with a precision of ±0.002 K in the 
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1.5 — 300 K range. The transition temperature Tq was 
determined from the maximum of the diffused intensity 
at small Q in unpolarized mode, with correction of the 
offset [14] due to the Fisher exponent r] and the finite Q 
value. The NSE spectra were recorded with a 0.45 nm 
neutron wavelength in the standard paramagnetic setup 
with a 2D multidetector whose pixels were divided in 
three constant-Q stripes. Separate measurements of the 
direct beam polarization made sure that any spurious de- 
polarization effects from the sample were avoided. The 
NSE spectra were measured for 1.001 Tc < T < 1.19 Tc 
and 0.43 < Q < 1.39 nm -1 , the scaling variable x = k/Q 
covering the range 0.03 — 1.6. The NSE resolution was 
measured on a similar spin glass sample at a sufficiently 
low temperature for the spin dynamics to be completely 
frozen and static on the probed timescale. All dynamic 
curves were divided by the corresponding resolution func- 
tion and normalized by the instantaneous magnetic re- 
sponse determined as part of the NSE measurements by 
x-y-z polarization analysis. This treatment gives the nor- 
malized intermediate scattering function Jjv(Qj^) that 
can be compared with the theoretical dynamical shape 
function Irg(Q, t) (equation [5]) . 
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FIG. 3: Dynamical scaling of the magnetic critical scattering 
at Tc in CdCr^S^ for Q > Qd- The red straight line displays 
the original RG model (equation [5]), the blue dash line a mod- 
ified RG model [15] and the green dash-dotted line a simple 
exponential decay of the magnetic correlations. 

At Tc, the adjustments of the Irg(Q,T) parametriza- 
tion to the NSE curves led to the power law decay (figure 
[TJ of the characteristic relaxation time r c = tqQ~ z with 
z = 2.47 and To = 0.02813 ns/nm~ z . Figure [3] shows 
that all spectra collapse in a master curve when plotted 
against the scaling variable t/r§Q~ z . The decay of the 
relaxation is faster than an exponential and is surpris- 
ingly well accounted for by the original RG calculation, 
even though it is only a first order e-expansion. A mod- 
ified RG model [15] is less convincing. We will comment 
on these modified RG versions later in the discussion. 

As mentioned previously, we went beyond the triv- 
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FIG. 4: Crossover of the shape function from the critical 
towards an exponential decay as a function of the reduced 
variable x = k/Q at x = 0.49 (upper figure) and at x = 0.77 
(the lower figure) . The blue solid lines display the RG theory 
at the corresponding x value, the red dash lines the critical 
shape and the green dash-dotted lines and exponential decay. 



ial verification of the dynamical scaling at Tc, where 
k/Q = 0, and measured the magnetic relaxation above 
Tc at finite values of x = k/Q, where the function de- 
parts from the critical shape and evolves towards an 
exponential. The spectra were measured at constant 
x = k/Q and scaled with the reduced time t/r(Q,x) = 
tQ(x)/roQ~ z . The experimental results are displayed in 
figure [4] for x = 0.48 and x = 0.77 respectively, with 
0(0.48) = 0.704 and 0(0.77) = 0.590. For each x, the 
magnetic relaxation spectra collapse on a single master 
curve, the shape of which evolves with increasing x to- 
wards the exponential. These master curves are well de- 
scribed by equation [5] with the measured critical expo- 
nents v = 0.7 ±0.03, z = 2.47 ±0.02 , and critical param- 
eters: transition temperature Tc = 84.84 ± 0.004 K and 
amplitudes ^0 = 3.61 nm -1 and To = 0.02813 ns/nm~ z . 

This exhaustive agreement between experiment and 
theory is new and in some sense surprising. Previous at- 
tempts to test the critical dynamic theory with neutron 
triple axis spectroscopy [10j[T5] adjusted the parameters 
a, k and b to the measured scaling function Q(x) over an 
extended x range, leading to large negative values for the 
parameter a. These induce a broad overdamped maxi- 
mum in F(Q,u) (figure [5| which is definitely not com- 
patible with experiments and our high resolution NSE 
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FIG. 5: Critical scattering F(Q, uS) of equation IT] for k = 0.51 
and —70 < a < —10, plotted for —0.5 < (JJIj c < 1 and 
showing an overdamped maximum for a < —10.75. 

experiments (figure [3]) when Fourier transformed in time. 
Obviously, these parameter adjustments are questionable 
as they are done in a much too large x range. In fact, 
only pure Iron follows the theoretical curve over such an 
extended range [8]. A self consistent reparametrisation 
would imply a one- step iteration to extract new uj c (Q, T) 
values with the reparametrized F(Q,u) and get a new 
scaling function Q(x) that would converge to a stable 
shape and set of parameters. Moreover, at large x values, 
only the b parameter might be adjusted [9] whereas the 
a parameter should be kept between 0.46 and 0.99 [16] 
or at least between -10.75 and 1 to avoid a broad maxi- 
mum in F(Q,u). The discrepancy between the neutron 
triple-axis results and the original RG predictions could 
therefore be used for a self consistent reparametrisation 
but even in this case the conclusions would be limited 
by the convolution of the complex function of equation [I] 
with the triple axis resolution function. 

Neutron Spin Echo is a unique technique for the de- 
tailed analysis of relaxation shape functions, especially 
magnetic ones. In NSE resolution corrections reduce to a 
simple division. Moreover, the paramagnetic NSE setup 
discards all nuclear parts of the scattering and extracts 
directly the critical intermediate functions, which are as 
exact as possible. This simple and straightforward data 
reduction allowed us to check the dynamical scaling on 
the shape of the relaxation over an extended temperature 
and Q range. Our NSE study shows that the RG theory 
gives an accurate quantitative description of the criti- 
cal scattering and dynamical scaling not only at x = 
but also at the crossover to the hydrodynamic regime at 
x < 1. At larger x values, however, the theory may not 
be applicable without reparametrization either due to its 
approximations or to deviations of the real systems from 
the 3D Heisenberg model. 

Similar studies were recently performed on Ising sys- 



tems showing an exponential decay on the microsecond or 
second timescales using Photon Correlation Spectroscopy 
with either X-rays [T7] or visible light [18]. However, the 
most detailed dynamic studies were performed by NSE 
[T9] on the nanosecond timescale. In spin glasses, dynam- 
ical scaling is more a question of phase space rather than 
reciprocal space, with streched exponential and power 
law decay [20] over an extended time range. The agree- 
ment between the experimental results and the theoreti- 
cal predictions presented in this letter goes beyond these 
findings. It is the first clear experimental validation of all 
the implications of the dynamical scaling hypothesis not 
only for the characteristic relaxation times but also for 
the very shape of the relaxation function, in the nanosec- 
ond and nanometer range. 
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